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Abstract 

Direct and indirect detection rates of relic neutralinos are reviewed in the frame- 
work of the Minimal Supersymmetric Standard Model. The theoretical estimates are 
compared with the most recent experimental limits from low-background detectors 
and neutrino telescopes. The properties of neutralino under the hypothesis that pre- 
liminary experimental results of the DAMA/Nal Collaboration may be indicative of 
a yearly modulation effect are examined. 
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ABSTRACT 

Direct and indirect detection rates of relic neutralinos are reviewed in the frame- 
work of the Minimal Supersymmctric Standard Model. The theoretical estimates 
are compared with the most recent experimental limits from low-background 
detectors and neutrino telescopes. The properties of neutralino under the hy- 
pothesis that preliminary experimental results of the DAMA/Nal Collaboration 
may be indicative of a yearly modulation effect are examined. 

1. Introduction 

Supersymmetric theories predict a large number of particles in excess to the Stan- 
dard Model ones. If the R-parity is conserved, the lightest among all the supersym- 
metric particles (LSP) must be stable. This feature makes the LSP a dark matter 
candidate, since this particle can be present today as a relic from the early stages of 
the evolution of the Universe. Different candidates have been proposed in the frame- 
work of supersymmetric theories: the neutralino or the sneutrino 1 in gravity medi- 
ated models, the gravitino 2 or some messenger fields in gauge mediated theories 3 , the 
axino 4 , stable non-topological solitons (Q-balls) 5 or others. 

In this paper we will focus on the most promising among all the different can- 
didates, the neutralino, which is defined as the lowest mass linear superposition of 
photino (7), zino (Z) and the two higgsino fields (H±, H%), i.e. x = a il + a 2^ + 
a^Hl + a^H^- The aim of this review is at providing the latest results on the cal- 
culation of different kinds of detection rates of relic neutralinos, in the framework of 
the Minimal Supersymmetric extension of the Standard Model (MSSM), constrained 
by the most recent experimental data coming from accelerator physics. We do not 
discuss here the details of the model, for which we refer to Refs. 6 ' 7 and to the ref- 
erences quoted therein. We only recall the standard assumptions employed here: i) 
all trilinear parameters are set to zero except those of the third family, which are 
unified to a common value A; ii) all squarks and sleptons soft-mass parameters are 
taken as degenerate: mj. = = mo; iii) the gaugino masses are assumed to unify 
at Mqut, and this implies Mi = (5/3) tan 2 Q^M 2 at the electroweak scale. After 
these conditions are applied, the free parameters are: M 2 , /i, tan/5, m^, m , A. The 
parameters are varied in the following ranges: 10 GeV < M 2 < 500 GeV, 10 GeV < 

^Report on the work done in collaboration with A. Bottino, F. Donato and S. Scopel. 
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Fig. 1. The scalar ncutralino-nucleon cross section ct^^ 01 ^, multiplied by the scaling factor £, 
is plotted versus the neutralino mass m x . The open curve denotes the 90% C.L. upper bound 
obtained from the total counting rates of Ref. 6 . The scatter plot represents the theoretical predictions 
calculated within the MSSM scheme (/U > only). The closed contour delimits the region singled 
out at 90% C.L. when the data of Ref. 12 are interpreted in terms of a modulation signal. 



< 500 GeV, 65 GeV < m A < 500 GeV, 100 GeV < m < 500 GeV, -3 < A < 
+3, 1.01 < tan/5 < 50. In our analysis the supersymmetric parameter space is con- 
strained by all the experimental limits on Higgs, neutralino, chargino and sfermion 
searches at accelerators. Moreover, the constraints due to the b — > s + 7 process 8 are 
satisfied. In addition to the experimental limits, we require that the neutralino is the 
lightest supersymmetric particle. Finally, the regions of the parameter space where 
the neutralino relic abundance exceeds the cosmological bound, i.e. Vt x h 2 > 1, are 
also excluded. 

2. Direct detection 

Neutralinos interact with matter both through coherent and spin dependent effects 
6,9 . We confine our discussion to the coherent effects, since these are the ones which 
are currently accessible to direct detection 9 . The relevant quantity, dependent on the 
supersymmetric parameters, which enters in the event rate of direct detection as well 
as in the indirect signals considered in Sect. 3, is p x x &i™ilr ° n \ where ^aiar °^ * s the 
X-nucleon scalar cross-section and p x is the neutralino local density. The expression 
of CTg°^ 1 r con ' ) and its relation to the differential event rate for elastic neutralino-nucleus 
scattering may be found in Ref. 9 . The x local density can be factorized as p x = £pi, 
i.e. in terms of the total local dark matter density pi. Here £ is evaluated as £ = 
min(l, fl x h 2 / (f2/i 2 ) m ; n ), where Q x h 2 is calculated as a function of the supersymmetric 
parameters 6 ' 10 and (Qh 2 ) min is a minimal value compatible with observational data 
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Fig. 2. Flux of up-going muons $^arth as a f unc ^i on f m ^ calculated within the MSSM scheme 
> 0). The solid (dashed) line represents the experimental 90% C.L. upper bound of Ref. 9 ( 10 ). 



and with large-scale structure calculations 7 . All the results of this paper refer to the 
choice pi = 0.5 GeV cm" 3 and {flh 2 ) min = 0.03 6 . 

Fig.l shows the present most stringent 90% C.L. upper limit on the quantity 
£ a lcailr 0n \ obtained by the DAMA/Nal Collaboration 11 . A possible rotation of the 
galactic halo can affect the exclusion plot less than a factor of two for m x > 50 GeV 12 . 
In the same figure our results are provided in the form of a scatter plot, obtained by 
varying the parameters of the supersymmetric model in the intervals quoted above 
(only fi > is displayed). We see that the scatter plot of the theoretical predictions 
reaches abundantly the curve of the 90% C.L. upper bound. This shows that the 
sensitivity of the direct detection experiment is adequate for a significant exploration 
of the neutralino parameter space. This result applies to the above mentioned choice 
of the astrophysical parameters pi and (f2/i 2 ) min . Different choices of these parameters 
inside their allowed ranges 6 ' 9 can change the theoretical predictions roughly by one 
order of magnitude. 

3. Indirect detection at neutrino telescopes 

Pair annihilations of neutralinos may provide other signals for the relic neutralino 
searches. The annihilation process may occur in the galactic halo or inside celestial 
bodies, like the Earth or the Sun, where the neutralino may have been accumulated 
as a consequence of gravitational capture. In the case of annihilation in the halo, the 
signal consists of photon 13 , positron and antiproton 14 fluxes, which can be observed 
by detectors placed on balloons or satellites. Recent calculations of the p fluxes in the 
MSSM will be presented elsewhere 15 . In the case of neutralinos captured in the Earth 
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Fig. 3. Flux of up-going muons <&^ un as a function of m x , calculated within the MSSM scheme 
(yu > 0). The solid (dashed) line represents the experimental 90% C.L. upper bound of Ref. 9 ( 10 ). 



and in the Sun, the signal is a flux of v^s which can be detected as up-going muons 
in a neutrino telescope. The values of the calculated fluxes of up-going muons 16 ' 17 
^Earth an( j ^Sun are displayed j n Figs.2,3 together with the experimental 90% C.L. 
upper bounds of Ref. 18 (solid line) and Ref. 19 (dashed line). From Figs.2,3 it turns 
out that many supersymmetric configurations give muon fluxes exceeding the present 
experimental upper bounds. Also in this case, the variation of the astrophysical 
parameters may change the calculated fluxes roughly by one order of magnitude. 

Another possible indirect signal, relevant to the neutrino telescope searches, may 
consist in a flux produced by neutrino annihilation inside the Large Magellanic 
Cloud (LMC). Theoretical estimates of the neutrino signal from LMC will be pre- 
sented elsewhere 20 . 

4. Direct detection: yearly modulation of the signal 

During the orbital motion of the Earth around the Sun, the direction of the 
velocity of the relic particles with respect to the detector chang function of 

time, and this induces a time dependence in the differential detection rate, i.e. 21 ' 22 
S(E,t) = So(E) + S m (E) cos[uj(t — to)], where ui = 27r/365 days and to = 153 days 
(June 2 nd ). Sq(E) is the average (unmodulated) differential rate and S m (E) is the 
modulation amplitude of the rate. The relative importance of S m (E) with respect 
to Sq(E) for a given detector, depends both on the mass of the dark matter particle 
and on the value of the recoil energy where the effect is looked at. Typical values 
of S m (E) j S (E) for a Nal detector are from a few percent up to roughly 15%, for 
neutralino masses of the order of 20-80 GeV and recoil energies below 8-10 KeV. 
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Fig. 4. The configurations compatible with the closed region of Fig. 1 are plotted in the m x -tan f3 
plane, within the gray area. Configurations which provide muon fluxes from the Earth and from 
the Sun which exceed the limits of Refs. 9,10 have been dropped. The dark region on the left side is 
excluded by current LEP data 3 . The region on the left of the vertical solid line will be accessible to 
LEP2 3 . The region on the left of the vertical dashed line will be explorable at TeV33 3 . In the region 
delimited by the closed dashed line, the neutralino relic abundance fl x h 2 may exceed the value 0.1. 



The DAMA/Nal Collaboration reported on an analysis of a collection of data over 
an exposure of 4549 Kg x days, obtained with an experimental set-up consisting of 
nine 9.70 Kg Nal(Tl) detectors 22 . The extraction of a possible signal is obtained by 
employing a maximum likelihood method applied to a binning in the recoil energy of 
the daily counts per detector. When interpreted in terms of a modulation signal due 
to a WIMP of mass m x and elastic cross section <y^^° n \ the data of Ref. 22 single out 
(at 90% C.L.) the region of Fig. 1 which is delimited by a closed contour (hereafter 
defined as region R m )- The occurrence of region R m as a domain relevant for a possi- 
ble modulation effect will necessarily require further investigation with much higher 
statistics. This point has clearly been stressed in Ref. 22 . Meanwhile, it is very inter- 
esting to analyze the possible implications that can be inferred from this preliminary 
result. Specifically, a number of questions deserve to be answered: a) what would 
be the features of a neutralino to satisfy the prerequisites of region R m ; b) would 
any other experimental search for relic neutralinos be able to investigate the region 
R m ] c) are neutralino configurations of region R m accessible to accelerator searches in 
the near future; d) what are the cosmological and astrophysical implications of relic 
neutralinos whose scalar elastic cross section is compatible with region R m 7 

We start our analysis by comparing the result of Ref. 22 with our calculation within 
the MSSM scheme. This comparison is shown in Fig.l, where we observe that many 
super symmetric configurations provide a value of £^sraiar° n ^ compatible with the re- 
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Fig. 5. The same configurations of Fig. 4 are displayed in the m/j-tan/3 plane, within the gray area. 
The dark regions are excluded by current LEP searches 3 or by theoretical arguments. Configurations 
which provide Q x h 2 > 0.1 fall within the region delimited by the closed dashed line. The region on 
the left of the solid line will be accessible to LEP2 3 . 



gion R m (we will call hereafter set S the supersymmetric configurations compatible 
with region R m ). As it was shown in Ref. 6 , many of these configurations provide 
sizeable muon fluxes from the Earth and the Sun, reachable by the foreseeable im- 
provements in neutrino telescopes. This result is interesting since it suggests the 
possibility to have an independent tool of exploration of the same configurations of 
set S. In the case of the indirect signal from the Earth, it is already possible to 
exclude part of the configurations of set S since they provide muon fluxes exceeding 
the Baksan and Macro Collaborations limits 6 . We therefore include in our following 
considerations only configurations which are not in conflict with indirect searches at 
neutrino telescopes (we define this subset as set T). 

Let us now consider the features of the supersymmetric configurations of set T. 
Fig. 4 shows, within the gray area, the neutralino mass ranges compatible with set T, 
for different values of tan f3. The dark area on the left side of the figure is excluded by 
current LEP data 8 . The region on the left of the vertical solid line around m x ~ 50 
GeV is the region explorable by LEP2 8 . We notice that, as far as the neutralino sector 
is concerned, LEP will be able to investigate only marginally the region of parameter 
space singled out by set T. A better chance to explore the neutralino mass range 
compatible with set T is given by the future upgrades at the Tevatron and by LHC. 
As an example, the region which extends up to the vertical dashed line is the mass 
region which will be possibly explored by TeV33 8 . 

Fig.5 shows the properties of the configurations of set T with respect to the lightest 
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Fig. 6. The neutralino local density p x , calculated for [^""^'^"^expt = 1 • 10 9 nbarn and pi = 0.5 
GeV cm~ 3 , is plotted versus the neutralino relic abundance fl x h 2 . For the value of [Scalar ° n ^expt 
employed here, the neutralino mass is restricted to the range 40 GeV < m x < 85 GeV, as obtained 
from the closed contour in Fig. 1. The two horizontal lines denote the physical range for the local 
density of non-baryonic dark matter. The two solid vertical lines denote the physical band for Qh 2 , 
and the two dashed lines give the preferred band for the cold dark matter contribution to tt. The 
two tilted dot-dashed lines denote the band where linear rescaling procedure for the local density is 
usually applied. 



Higgs mass m h , a parameter which is crucial in establishing the size of both direct 
and indirect detection signals. The gray area in Fig. 5 corresponds to values of rrih 
and tan (3 compatible with set T. The dark regions are excluded by current LEP 
searches 8 (on the left of the plot) or by theoretical arguments (on the right side). In 
Fig. 5 it is also reported the region which will be accessible to LEP2 8 . A large portion 
of the region compatible with the modulation analysis will be covered by the LEP2 
searches. In particular, all the region for tan j3 < 3 will be analyzed. 

In Ref. 6 it was shown that a large number of configurations belonging to set T 
provide a large value of the neutralino relic abundance. These configurations are 
very appealing from the point of view of dark matter, and they are shown in Fig.4,5 
where the closed dashed line contains the portion of the parameter space where Q x h 2 
may exceed the value 0.1. In order to investigate the cosmological and astrophysical 
properties of the configurations of set T, we present an analysis which is meant 
to obtain from the experimental data the relevant cosmological implications for relic 
neutralinos in the most direct way. We adopt the following procedure 7 : 1) we evaluate 
a scaiar C ° n ' ) an d ^x^ 2 varying the supersymmetric parameters; 2) for any value of 
[Px°"s ( caiar° n) ]cxpt = Pi [C^scai^^lexpt compatible with the experimental region R m we 
calculate p x = [Px^scaiar^lexpt/ ^scalar °^ an d restrict the values of m x to stay inside 



the region R m . 3) The results are displayed in a scatter plot in the plane p x vs. Q x h 2 . 

One example of our results is given in Fig. 6 for [C^icaiax^lexpt = 1 • 10~ 9 nbarn 
(pi = 0.5 GeV cm~ 3 is used). The two horizontal lines denote the physical range 0.1 
GeV cm~ 3 < p\ < 0.7 GeV cm" 3 for the local density of non-baryonic dark matter 7 . 
The solid vertical lines denote the physical band for Qh 2 : 0.03 < Qh 2 < 1, and the 
two dashed lines give the favored band for the cold dark matter contribution to Q: 
0.1 < (fi/i 2 ) C DM < 0.3 7 . The two tilted dot-dashed lines denote the band where 
linear rescaling procedure for the local density is usually applied. With the aid of this 
kind of plot we can classify the supersymmetric configurations belonging to region 
R m into various categories. Configurations whose representative points fall above the 
maximum value p x = 0.7 GeV cm -3 have to be excluded (those providing an Q x h 2 > 1 
are already disregarded in the plot). Among the allowed configurations, those falling 
in the region inside both the horizontal and solid vertical lines are very appealing, 
since they would represent situations where the neutralino could have the role of a 
dominant cold dark matter component; even more so, if the representative points fall 
in the subregion inside the vertical band delimited by dashed lines. Configurations 
which fall inside the band delimited by the tilted dot-dashed lines denote situations 
where the neutralino can only provide a fraction of the cold dark matter both at the 
level of local density and at the level of the average f2. Configurations above the upper 
dot-dashed line and below the upper solid horizontal line would imply an unlikely 
special clustering of neutralinos in our halo as compared to their average distribution 
in the Universe. 

5. Conclusions 

In this paper we have reported the most recent calculations of the direct and in- 
direct detection rates of relic neutralinos in the framework of the Minimal Supersym- 
metric Standard Model. We have shown that the theoretical estimates of the detection 
rates may be at the level of the present experimental sensitivities of low-background 
detectors and neutrino telescopes. For many supersymmetric configurations, and for 
median values of the astrophysical parameters which enter in the calculations of the 
detection rates, the predicted signals may already exceed the present experimental 
bounds. This shows that the different experimental efforts to search for relic particles 
are potentially able to deeply investigate the possibility that neutralino is a compo- 
nent of the dark matter of the Universe. An interesting preliminary analysis of the 
DAMA/Nal Collaboration has shown that their data are compatible, at 90% C.L., 
with a modulation signal of the direct detection rate. The features of a neutralino 
able to satisfy the prerequisites of this signal have been analyzed and it has been 
shown that many configurations are compatible with a dark matter scenario where 
the neutralino is the major component, both on galactic and cosmological scales. 
However, we have to remind here that the occurrence of a possible modulation ef- 
fect will necessarily require further investigations with much higher statistics. This 



project is currently under way. 
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